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Abstract. The Gilbert type bound for codes in the title is reviewed, both 
for small and large alphabets. Constructive lower bounds better than these 
existential bounds are derived from geometric codes, either over Fp or Fp2, or 
over even degree extensions of Fp. In the latter case the approach is concatena- 
tion with a good code for the Hamming metric as outer code and a short code 
for the Lee metric as an inner code. In the former case lower bounds on the 
minimum Lee distance are derived by algebraic geometric arguments inspired 
by results of Wu, Kuijper, Udaya (2007). 



1. Introduction 

The Lee metric was introduced in the coding theoretic hterature as a way to 
analyze codes over large alphabets applied to phase modulation More recently, 
it was used successfully in the correction of errors in constrained memories [7] . The 
main challenge in Lee codes has always been to generalize the known results for the 
Hamming metric at the price of added technical complexity. For instance the ana- 
logues of the classical bounds of Hamming and Gilbert that are derived in [5] and 
are reviewed below do not even admit a closed form. In the present work we will 
strive to generalize what is perhaps the most unexpected and most difficult result 
of algebraic coding theory in the last century: the Tsfasman-Vladut-Zink bound [5] 
that shows the existence of families of codes strictly better than the Varshamov- 
Gilbert bound, a non-constructive existential bound. 

We begin by reviewing the known analogues of the Gilbert bound in the Lee 
metric (we reserve the term Varshamov bound for algorithmic constructions of gen- 
erator matrices). The analogue of the Gilbert bound was derived first for odd 
alphabet size by Astola [2], using Lagrange multipliers, and for all q by Gardy- 
Sole |3] using saddle point approximation. In the same paper a bound valid for 
so-called large alphabets (say q > 2n + 1) was introduced. All these bounds are 
non-constructive; in the present note we compare them with constructive bounds 
obtained from algebraic geometric constructions. In particular a bound on the Lee 
minimum distance of geometric codes over a prime field was derived recently in [9] . 
Unfortunately the authors of this interesting paper did not take into account the 
case of constant functions on the curve, and did not consider asymptotics. We cor- 
rect this small technical oversight and complete their work by deriving asymptotic 
estimates of the parameters of the shortened codes. We also derive a variant of 
this bound over Fp2 that allows to use geometric codes on the TVZ bound when 
the Ihara function A{q) is known exactly A(p'^) = p — 1. In general estimating 
A{p) for p a prime is a difficult question [6j. This latter bound is compared with 
the natural idea of concatenating a good Hamming weight geometric code over an 
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extension field with a BCH code over the base field with a good Lee distance [7]. 
Plotting graphs for asymptotic bounds show that the corrected bound of [3] is in 
general better than the concatenated code bound, which in turn outperforms the 
non-constructive Gilbert type bound. 



2. Definitions 

If q is any natural number we recall that the Lee weight of a symbol x of Zg is 
min(x, q — x). The Lee weight of a vector is then defined by adding the contribution 
of each symbol. 

The geometric codes we will consider will be of the form C{D,G) with D,G 
disjoint divisors of an algebraic curve of genus g over F^. Here D = Pi + • • • + P„, 
is a localizing set; a possible choice for G is rP, for some P on the curve; and the 
codeword attached to a function / in the L{G) space is of the form 

(/(Pi),...,/(F„))- 

We let r = deg(G), and assume throughout that r > 2g — 2 so that the parameters 
of C{D, G) are [n, k] with k = r — g + 1. By the Ihara function A{q) we shall mean 
the largest number of rational points of a curve over Fg per unit of genus, for genus 
going to cx). In symbols 

A(g)-limsup^^, 

g-y+oo g 

where Nq{g) is the largest number of rational points of a curve over Fg, also known 
as the Serre function. It is known that A{q) < ^ — 1, for all q (Drinfeld/Vladut 
bound) and by using modular curves or recursive towers that A(p^"^) = p"^ — 1. 
The coding motivation to define the Ihara function is the so called TVZ bound. 
This asymptotic bound states 

R+S>1 ^ 



A{qy 

for R the rate of the family of geometric codes considered and S its relative distance. 

3. Gilbert bound 

Let R{6) be the largest achievable rate of a family of codes of relative Lee distance 

6. 

3.1. Small alphabets. The following bound was obtained in [21 Theorem 2]. 
Theorem 3.1. If q = 2s + 1, then R{6) > 1 + log^ where a, /3 are defined by 
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3.2. Large alphabets. The following function plays the role of Shannon entropy 
in the present context. 



Lg{x) = x\ogq X + logq{x + Va;2 + 1) - x\ogq{y^x'^ + 1 - 1). 
The following result was first derived in [3] . 
Theorem 3.2. If q > 2dn + 1 and S < L-\l) w 0.37, then 

R{6) > 1 - L{6). 

4. Geometric codes 

4.1. Concatenation. By the TVZ bound [Si Theorem 13.5.4] we know there are 
families of geometric codes over GF{Q) for Q a square, with rate TZ and relative 
distance A satisfying 

7^ + A > 1 - —1 

We concatenate this geometric code with a code over Zq of parameters [n, k] and 
minimum Lee distance d^. We must assume therefore that Q = . If g is not a 
prime, we label Zg by the elements of GF{q) in an arbitrary fashion. In order to 
apply the TVZ bound we must assume Q to be a square, or, equivalently k to be 
even. 

Proposition 4.1. With the above notation the rate R and the relative Lee distance 
S of the concatenated code satisfy 

^ ' >l 1 



k/n d^/n q^/'^ — 1 

Using some BCH codes over Zp as inner codes yields the following bound. 

Corollary 4.2. For each prime p > 7 and every integer 1 < t < (p + l)/2, such 
that p is congruent tot+1 mod 2, there is a family of Lee codes overZp with rate 
R and relative Lee distance S satisfying 

Rjp^l) S{p-1) 1 
p-l-t 2t - p(p-*-i)/2 - 1' 



Proof. For these inner codes n ^ p — 1, d^ > 2t and /c > n — i, by [7]. □ 

4.2. Victorian bound. In [9j a Lee analogue of Goppa estimate was derived by 
three authors from Victoria state of Australia. There was a slight error in their 
bound but it can easily be fixed as follows: 

Theorem 4.3. Given an algebraic curve of genus g over ¥q having at least n + 1 
rational points, there are codes of parameters [n~l,r — g] overWq with Lee distance 

dj, > 



Ar 

for any integer r in the range {2g — 2, n). 
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Proof. The geometric codes C{D,rP) of j9] have non-constant codewords of Lee 

2 _ 2 

weight > " , and constant codewords of Lee weight n, 2n, ... To remove these 
constant codewords we shorten the code in an arbitrary position thus decreasing 
the dimension by one unit. □ 

The following Corollary is immediate. 

Corollary 4.4. For a family of curves of genus g ^ 771, the rate R of the attached 
family of codes of relative distance S is 

i? > -7 - 25 + v/4(52 + 1. 



Proof. Assume that for n — > cx), we have r xn; then passing to the limit in 
k = r — g, we get a; = i? + 7. Using the above theorem we should have 

P{x) = + 45a; - 1 > 0, 

which happens for x > only if x is larger than the larger root of P{x) = 0. □ 

The best choice of 7 is 7 = 1/A{q)^ in particular if g is a square one can take 
7 = ~ 1) as in the TVZ situation, but the resulting code may be non- linear 

over the ring Zg. If g = p is a prime, less is known about the A{p)^ but for example 
from we know that A(5) > 0.727. 

4.3. A new bound using descent of the base field. Let p be an odd prime. 

Let {1, a] be a basis of Fp2 over Fp, so Fp2 = Fp-l + Fp-a = Fp xFp. If c g (Fp2)" is 
a word of length n over Fp2 , then using this identification we get a word c e (Fp)^" 
of length 2n, and likewise, a linear code C of parameters [n, k] over Fp2 gives rise to 
a linear code C of parameters [2n, 2k] over Fp. We extend the definition of the Lee 
weight to Fp2 by setting the weight of a symbol z = x + ya ^ Fp2 (where x,y Fp) 
as 

wtL^z) = wtL(x) + wtL{y), 
so summing over coordinates we find 

wtL(c) = wtL{c). 

We also let S{j) be the cardinality of the "sphere" of radius j in Fp2 , so 

S{j)^\{ze¥p2; wtL{z)^j}\ 
and B{M) the cardinality of the "ball" of radius Af, 

M 

B{M) = \{z e Fp2; wtLiz) < M}\ = ^^(j). 

Last, we let W{t) be the sum of the weights of the t smallest elements in Fp2, so 
formally 

M 

= ^jS'(j) + (M + l)(t-B(M)) forB(M) <t<B{M + l). 
3=0 

In this formula we allow real values of t. 

These quantities can be computed explicitly. For 1 < j < the symbols 
z S Fp2 having Lee weight j are precisely the z = ±i ± {j — i)a for < i < J, and 
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Figure 1 . Lee spheres and Lee balls for F72 = F7 x F7 
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W is linear between W{B{M) and W{B(M + 1)) 

Figure 2. Values of the functions S, B, and W, for p = 7 

for < j < p—1 these are the z — ±{^^ —i)±{j ~ ^^+i)a for < i < p— 1 — j, 
so 



1 



S{j) = { 4j 



if J = 

if 1 < J < ^ 



4{p-j) ifH±i<j<p-l 



from which it follows 



5(M) JV'''^''^'^ 



if < M < 



p-i 



- 2{p-M){p-l-M) if 2±i < M < p - 1 
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and then, for B{M) <t< B{M + 1), a straightforward computatfon gives 

W{t) = 

( 2M(Af+l)(2M+l) ^ ^M + l){t-B{M)) if M < 

j p(p-l)(p+l) „ 2(p-M)(p-l-M)(p+l+2M) ^ (M + l)(t-i3(Af)) ifAf>2+i 

which can also be written 

j (Af+l)(i-l) - if M < 2=i 

\(Af+l)ft-p^) + - 2(p-M)(p^l-M)(p-2-M) if^,^>P±l, 

For example, we see that W has maximal value 

Also, remark that in the case M < the expression does not depend on p. And 
in the particular case t = B{M) for M < 2=^ we have t = 1 + 2M{M + 1) so 
M = i(^/2^^ - 1), so then W{t) = |Af (A-/ + l)(2Af + 1) = M^approa;(i) where 

1 

3* 



Wapproxit) = Tj(i - 1)V2<-1. 



* * >^ »'((),p = 3 



X X X m),p'5 

^ WW.p = 7 




Figure 3. Graphs of Wapprox{t) and W{t) for p = 3, 5, 7 

In fact one can show that Wapprox is a lower bound for W: 
Lemma 4.5. We have 

W{t) > Wapproxit) 

for all real t > 1. 
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Proof. First define Bo(Af ) = l+2il/(A/+l) for all integer M > 1, and then Wo{t) = 
|M(M + l)(2A-f + 1) + (M + l)(t-l-2M(M+l)) for Ba{M) <t< Bo{M + l). 
By construction we have B{M) = Bq{M) for M < and B{M) < Bo{M) 

for M > 2zJ:. From this it follows W{t) = Wo(i) for I < t < p^+^p-^ ^ and 

W{t)>Woit)ioTt> P'+f-\ 

Computing its second derivative, we see Wapprox{t) is a convex function of < > 
1/2. On the other hand Wo(t) is a piecewise linear function, and by construction 
both take the same values for all t of the form t = 1 + 2M(M + 1). Since a convex 
function lies below its chords, we then conclude W{t) > Wo{t) > Wapprox{t) for all 
real t>l. □ 

It is interesting to note that since Wapproxit) is a convex interpolation of the 
piecewise linear function W{t) in the range 1 < t < ^ ^ , not only will it be a lower 
bound, but also a reasonably good approximation for it. However, note that out of 
this range, while it is still a lower bound, the approximation can be slightly worse. 
We compare the graphs of Wapproxit) and W{t) for some small primes in Figure |3] 

Theorem 4.6. Consider an algebraic curve of genus g over Fp2 having at least 
n + 1 rational points, and an integer r in the range {2g — 2,n). Let M be an integer 
such that B{M) < n/r < B{M + 1), e.g. 



M 



i(-l + v/2n/r-l)J ifl< n/r < p"+f-3 

+ v/V + l-2nA)J P'+f-^ < n/r < p\ 



Then there are codes of parameters [2(n — 1), 2(r — g)] over the prime field ¥p with 
Lee distance 
dL > rW{n/r) 

~ j(M + l)n+ r */n/r>H^±|H^ 



> 



rWapprox{n/r) = 



2n 



Proof. As before let D = Pi + • • • + P„ and G = rP where P,Pi, . . . ,Pn are rational 
points on the curve, so that C = C{D, G) has parameters [n, r + 1 — g] over Fp2 , 
and C has parameters [2n, 2(r + 1 — g)] over Fp. 

Now following the idea from [9] (already present in [7]), we remark that if c is a 
non-constant word in C, then any symbol from Fp2 can occur at most r times in c. 
This implies that we will have wt]^{c) > wt^ia) where a = (ai, 02, . . . , a„) G (Fp2)" 
is constructed as follows: 

• all symbols z e Fp2 of Lee weight 0, 1, . . . , M occur in a exactly r times 
each 

• some symbols of Lee weight M -\- 1 could occur in a, but not more than r 
times each, and at least one of them less than r times 

• no symbol of Lee weight greater than M + 1 occur in a 

where M is such that B{AI) < n/r < B{M + 1). By construction we then have 

wtL^a) — rW(n/r). 
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To conclude we shorten C to get rid of the constant codewords, decreasing the 
length and the dimension of C by one, and the length and the dimension of C 
by two, so this new C has parameters [2n — 2, 2(r — g)] over Fp and Lee distance 
(11 > rW{n/r), as claimed. □ 

Corollary 4.7. Suppose we are given a family of curves over Fp2 having at least 
n + 1 rational points and genus g ^ ^n, for some integers n — >■ oo. Then from 
this, for any R < I — j, we get a family of codes over Fp of asymptotic rate R and 
asymptotic relative Lee distance 



6 V ^ + 7 



^ > < 



Proof. Apply Theorem 4.6 with r ^ xn where x — R + j and pass to the limit. □ 

Remark we have 7 > hence in the proof n/r ^ 1/x < p — 1 < ^ < 

P-+^Li!L^ so^ of the two possible expressions for W we only need the one corre- 
sponding to this range. As already noted, this expression does not depend on p, 
and moreover, this is the range where Wapprox is a good approximation of W . 
The first inequality in this corollary can then be put in the form 

fi(l-(i? + 7)) iil<T^<5 {p>3) 

i(2-6(i? + 7)) if5<^<13 (p>5) 

i(3-19(i? + 7)) ifl3<^<25 {p>7) 
etc. 

, AM - bu {R + 7) if BiM) <j^<BiM + l) {p> 2M + 3) 

where aj\/ = ^'^^^ and = (^-f+i)(2M +4M+3} ^ second inequality is a 

reasonably good approximation of the first and might be easier to handle in some 
applications, especially when i? + 7 is small. 

We conclude with an equivalent formulation of the previous corollary. Denote 

by 

aLee,p{5) 

the largest real R such that there exists a family of codes over Fp of asymptotic 
rate R and asymptotic relative Lee distance S. Introduce the function 

2 X 

and 

X 1 
g{x) = -Wapproxi-)- 

Then: 

Corollary 4.8. Letj^^. We have the lower bounds 



lower bounds on the minimum distance of long codes in the lee metric9 

Table 1. 





5^ 




11^ 


13^ 


17^ 


19^ 


23^ 




0.0308 


0.0095 


0.0023 


0.0014 


0.0007 


0.0005 


0.0003 



These inverse functions / ^ and g ^ can be computed explicitly. This gives the 
expressions: 



1-26 -J 
i(l-J)-7 

2S) 



if < J < I (p > 3) 
if 1 < (5 < i§ (p > 5) 
(P>7) 



5 

10 



28 



19 V-^ 3 J I 13 — " — 25 

etc. 

CM - dMS - 7 if C{M) <S< C{M + 1) {p> 2M + 3) 



where cm = 2A-P+4M+3 ' "a^ 
and 



(A/+l)(2M2+4J\/+3) ' 



C(M) 



— M+1 (A/+l)(2M^+4M+3) 
2 6(l+2Af(A//+l)) ' 



f|cos ^( 



1/3 



if A > 



2i/ cos I 



>7v'^ \ I 27r 
-4m3 ) 3 



I - 7 if A < 



with A = 6912(56 _^ 2112(5* - i^, w = 365^ - i and w = (48(5^ - |) 



5. Comparisons 

5.1. Concatenation vs Victoria and descent. For large p and t of order p/2, 
the bound of Corollary |4.2| can be approximated by 

2i? + (5 = 1, 

For large q, i.e 7 = 0, the bound of Corollary |4.4| can be approximated by 



R>-25+ V4(52 + 1. 
Clearly the latter is above the former. 

For finite values of p in Corollary |4.2| and q the square of a prime in Corollary |4.4[ 
this first happens for p — 113 or 127 and q — 11^ (sa p) as shown in Figure [4j 
For q — p being prime, Corollary |4.2| is improved by Corollary 4.7 where the case 
for (7 = 7 is as shown in Figure [6j 



5.2. Astola vs Victoria and descent. For q at least 5^ Corollary 



4.4 



IS an im- 



provement of Theorem 3.1 for S > 5q ~ 0.0308. We list some small values of q and 
such Sq in Table [l] as well as give a graph of comparison for q = 23^ in Figure 5] 
For q being prime, Theorem |3.1| is improved by Corollary |4. 7| where the case for the 
first prime g = 1 1 is as shown in Figure [7] 
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Concatenation (p^113) 
Concatenation (p=127) 
Victoria (q=121) 
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Figure 4. 
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Figure 7. 



